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Abstract 
Bau = 0 B29 > 0 Boo < 0 Bao > 0 
In this contribution we present some results of potential energy sur- 
faces of actinide and transfermium nuclei from multi-dimensional 
= constrained relativistic mean field (MDC-RMF) models. Recently 
we developed multi-dimensional constrained covariant density func- 
* tional theories (MDC-CDFT) in which all shape degrees of freedom Boo #0 B30 #0 B32 #0 Boo > 0 


~ Êu with even p are allowed and the functional can be one of the 
following four forms: the meson exchange or point-coupling nucleon 
a interactions combined with the non-linear or density-dependent cou- 


) plings. In MDC-RMF models, the pairing correlations are treated 


“= with the BCS method. With MDC-RMF models, the potential 


energy surfaces of even-even actinide nuclei were investigated and 
the effect of triaxiality on the fission barriers in these nuclei was 
= discussed. The non-axial reflection-asymmetric 632 shape in some 
transfermium nuclei with N = 150, namely 246Cm, 248Cf, 25°Fm, 


and 252No were also studied. 


1. Introduction 


The occurrence of spontaneous symmetry breaking in 
) atomic nuclei leads to various nuclear shapes which can 
usually be described by the parametrization of the nu- 
clear surface or the nucleon density distribution [1,2]. In 
mean-field calculations, the following parametrization 


AT 
~ ZARA (1) 


is usually used, where @), are the mass multipole op- 
erators. In Fig. 1, a schematic show is given for some 
typical nuclear shapes. The majority of observed nuclear 
shapes is of spheroidal form which is usually described by 
B20. Higher-order deformations with A > 2 such as (39 
also appear in certain atomic mass regions [3]. In addi- 
tion, non-axial shapes in atomic nuclei, in particular, the 
nonaxial-quadrupole (triaxial) deformation (622 has been 
studied both experimentally and theoretically [4-6]. The 
influence of the nonaxial octupole (32 deformation on the 
low-lying spectra has been also investigated [7-16]. 
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Fig. 1: (Color online) A schematic show of some typical nuclear 
shapes. From left to right, the 1st row: (a) Sphere, (b) Prolate 
spheroid, (c) Oblate spheroid, (d) Hexadecapole shape, and 
the second row: (e) Triaxial ellipsoid, (f) Reflection symmetric 
octupole shape, (g) Tetrahedron, (h) Reflection asymmetric 
octupole shape with very large quadrupole deformation and 
large hexadecapole deformation. Taken from Ref. [17]. 


In nuclear fission study, various shape degrees of free- 
dom play important and different roles in the occurrence 
and in determining the heights of the inner and outer bar- 
riers in actinide nuclei (in these nuclei double-humped fis- 
sion barriers usually appear). For example, the inner fis- 
sion barrier is usually lowered when the triaxial deforma- 
tion is allowed, while for the outer barrier the reflection 
asymmetric (RA) shape is favored [18-22, 24, 25]. 

In order to give a microscopic and self-consistent de- 
scription of the potential energy surface (PES) with more 
shape degrees of freedom included, multi-dimensional con- 
strained covariant density functional theories were devel- 
oped recently [26,27]. In these theories, all shape degrees 
of freedom àu with even ys are allowed. In this contribu- 
tion, we present two recent applications of these theories: 
the PES’s of actinide nuclei and the non-axial reflection- 
asymmetric (632 shape in some transfermium nuclei. In 
Section 2, the formalism of our multi-dimensional con- 
strained covariant density functional theories will be given 


briefly. The results and discussions are presented in Sec- 
tion 3. Finally we give a summary in Section 4. 


2. Formalism 


The details of the formalism for covariant density func- 
tional theories can be found in Refs. [28-33]. The CDFT 
functional in our multi-dimensional constrained calcula- 
tions can be one of the following four forms: the meson 
exchange or point-coupling nucleon interactions combined 
with the non-linear or density-dependent couplings [26, 
27, 34, 35]. Here we show briefly the one corresponding to 
the non-linear point coupling (NL-PC) interactions. The 
starting point of the relativistic NL-PC density functional 
is the following Lagrangian: 


L= (iy ,o" Myy Liin Lol Laer Loou; (2) 
where 
1 1 1 1 
Lin = agp + ~avpy + =arsprs + -arv pry, 
2 2 2 2 
1 1 1 
La = zôsPs + zasos Er lev]? 
1 1 1 23 
Laer = 55s |Oves])” + 30v [3v pv]? + zôrslðv prs] 
1 - 
+ 3 oTV [3v örv]? 
1 p: 1— T3 
Loou = re Pup + e— Aopv, (3) 


AJ are the linear coupling, nonlinear coupling, derivative cou- 
‘= = pling, and the Coulomb part, respectively. M is the nu- 
=> cleon mass, as, ay, ars, ATV, Bs, Ys, YV, Os, Ôv, Ors, 
àg and dry are coupling constants for different channels and 
em ¢ is the electric charge. ps, prs, pv, and pry are the 
¢ = isoscalar density, isovector density, time-like components 
= of isoscalar current, and time-like components of isovec- 


tor current, respectively. The densities and currents are 


© defined as 


prs = Ty, 
Prv = YPY. 


PS = pp, 
pv = 07, (4) 
Starting from the above Lagrangian, using the Slater de- 
terminants as trial wave functions and neglecting the Fock 
term as well as the contributions to the densities and cur- 
rents from the negative energy levels, one can derive the 
equations of motion for the nucleons, 


hab = (a -P+ B(M + S(F))+V(F)) vi = edi, (5) 


where the potentials V(r) and S(r) are calculated as 


S = agsps+Bspst+ ses + 6sAps 
+ (arsprs + 6rsAprs) T3, (6) 
V = ayvpy + wey + dvApvW 


+(arvprv + orvAprv) T3. 


An axially deformed harmonic oscillator (ADHO) ba- 
sis is adopted for solving the Dirac equation [26, 27, 36]. 
The ADHO basis is defined as the eigen solutions of the 
Schrodinger equation with an ADHO potential [37, 38], 


R2 
—s7V? + Ve(2,p)| ®alro) = 


where 1 
Va(z, p) = 5M wpe a aden", 


(9) 


is the axially deformed HO potential and w; and wọ are the 
oscillator frequencies along and perpendicular to z axis, 
respectively. These basis states are also eigen functions 
of the z component of the angular momentum j, with 
eigen values K = mı + ms. For any basis state ®a(ro), 
the time reversal state is defined as Palro) = T®,(ro), 
where T = io, is the time reversal operator and K is 
the complex conjugation. Apparently we have Kg = -Ka 
and ta = Ta. These basis states form a complete set for 
expanding any two-component spinors. For a Dirac spinor 
2a fi Pa(re) 


with four components, 
Va Palro) ) 


wi(ro) = ( 

where the sum runs over all the possible combination of 
the quantum numbers a = {nz,n-,mi,ms}, and f7 and 
g& are the expansion coefficients. In practical calculations, 
one should truncate the basis in a proper way [26, 27, 36]. 

The nucleus is assumed to be symmetric under the V4 
group, that is, all the potentials and densities can be ex- 
panded as 


(10) 


1 = 1 
F(p, p2) = folo 2) 7 + d, Sales) Fe cos(2ny), 
(11) 


The PES is obtained by the constrained self-consistent 
calculation, 
1 
E' = Epur + ` OQA» (12) 
Au 


where the variables C),,’s change their values during the 
iteration. 

Both the BCS approach and the Bogoliubov transfor- 
mation are implemented in our model to take into ac- 
count the pairing effects. For convenience, we name the 
MDC-CDFT with the BCS approach for the pairing as the 
MDC-RMF model and that with the Bogoliubov trans- 
formation as the MDC-RHB model. More details of the 
multi-dimensional constraint covariant density functional 
theories can be found in Refs. [26, 27]. 


3. Results and discussions 


3.1. PES’s of actinides 


In Refs. [26,27], one- (1-d), two- (2-d), and three- 
dimensional (3-d) constrained calculations were performed 
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Fig. 2: (Color online) Potential energy curves of 74°Pu with 
various self-consistent symmetries imposed. The solid black 
curve represents the calculated fission path with V4 symme- 
try imposed: the red dashed curve is that with axial symmetry 
(AS) imposed, the green dotted curve that with reflection sym- 
metry (RS) imposed, the violet dot-dashed line that with both 
symmetries (AS & RS) imposed. The empirical inner (outer) 
barrier height Bemp is denoted by the grey square (circle). The 
energy is normalized with respect to the binding energy of the 
ground state. The parameter set used is PC-PK1. Taken from 
Ref. [26]. 


for the actinide nucleus 74°Pu. The MDC-RMF model 
with the parameter set PC-PK1 [39] was used. In Fig. 2 we 
show the 1-d potential energy curves from an oblate shape 
with 629 about —0.2 to the fission configuration with 820 
beyond 2.0 which are obtained from calculations with dif- 
ferent self-consistent symmetries imposed: the axial (AS) 
or triaxial (TS) symmetries combined with reflection sym- 
metric (RS) or asymmetric cases. The importance of the 
triaxial deformation on the inner barrier and that of the 
octupole deformation on the outer barrier are clearly seen: 
The triaxial deformation reduces the inner barrier height 


' by more than 2 MeV and results in a better agreement 


with the empirical value [22]; the RA shape is favored be- 
yond the fission isomer and lowers very much the outer 
fission barrier [21]. Besides these features, it was found 
for the first time that the outer barrier is also consider- 
ably lowered by about 1 MeV when the triaxial deforma- 
tion is allowed. In addition, a better reproduction of the 
empirical barrier height can be seen for the outer barrier. 
It has been stressed that this feature can only be found 
when the axial and reflection symmetries are simultane- 
ously broken [26]. 

Two-dimensional PES’s in the (6209, G22) plane near the 
inner and outer barriers are shown in Figs. 3 and 4, re- 
spectively. Starting from the axially symmetric ground 
state, the nucleus goes through the triaxial valley to the 
isometric state. The inner barrier is located at 629 ~ 0.65 
and {$22 ~ 0.06. The isomeric state keeps an axially sym- 
metric shape. As (29 further increases, the nucleus goes 
through a triaxial valley again, and then goes fission. The 
outer barrier is located at Boo œ% 1.21, Boo œ% 0.02, and 
B30 & 0.37. 


Fig. 3: (Color online) Potential energy surfaces of 7*°Pu in 
the (G20, G22) plane around the inner barrier. The energy is 
normalized with respect to the binding energy of the ground 
state. The least-energy fission path is represented by a dash- 
dotted line. The saddle point is denoted by the full star. The 
contour interval is 0.5 MeV. 
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Fig. 4: (Color online) Potential energy surfaces of 74°Pu in 
the (G20, G22) plane around the outer barrier. The energy is 
normalized with respect to the binding energy of the ground 
state. The least-energy fission path is represented by a dash- 
dotted line. The saddle point is denoted by the full star. The 
contour interval is 0.25 MeV. 


EN) proton numbers [9, 12, 16]. 
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Fig. 5: (Color online) The binding energy E (relative to the 
ground state) for °48Cf as a function of the non-axial octupole 
deformation parameter (32. 


A systematic study of even-even actinide nuclei has been 
carried out and the results were presented in Ref. [27] 
) where we have shown that the triaxial deformation lowers 
) the outer barriers of these actinide nuclei by about 0.5 ~ 1 
) MeV (about 10 ~ 20% of the barrier height). 


3.2. Y32-correlations in N = 150 isotones 


It has been anticipated that the tetrahedral shape (Bry = 
0, if A Æ 3 and u Æ 2) may appear in the ground states of 
) some nuclei with special combinations of the neutron and 
The tetrahedral symmetry- 
* driven quantum effects may also lead to a large increase 
of binding energy in superheavy nuclei [40]. However, no 
solid experimental evidence has been found for nuclei with 


© tetrahedral shapes. On the other hand, (32 deformation 


may appear together with other shape degrees of free- 
dome, say, (2. For example, it has been proposed that 
the non-axial octupole Y32-correlation results in the exper- 
imentally observed low-energy 27 bands in the N = 150 
isotones [41] and the RASM calculations reproduces well 
the experimental observables of these 27 bands [42]. 

In Ref. [43] the non-axial reflection-asymmetric $32 de- 
formation in N = 150 isotones, namely 74°Cm, 745Cf, 
250Fm, and 7°2No was investigated using the MDC-RMF 
model with the parameter set DD-PC1 [44]. It was found 
that due to the interaction between a pair of neutron or- 
bitals, [734]9/2 originating from vj15/2 and [622]5/2 orig- 
inating from vgg/2, and that of a pair of proton orbitals, 
[521]3/2 originating from mf7/2 and [633]7/2 originating 
from 7713/2, rather strong non-axial octupole Y39 effects 
appear in 748Cf and 2°°Fm which are both well deformed 
with large axial-quadrupole deformations, 820 ~% 0.3. 

In Fig. 5, the potential energy curve, i.e., the total bind- 
ing energy as a function of 832 was shown for 248Cf. At 
each point of the potential energy curve, the energy is min- 
imized automatically with respect to other shape degrees 
of freedom such as (20, 222, 630, and 840, etc. One finds 


in this curve a clear pocket with the depth more than 0.3 
MeV. Similar potential energy curve was also obtained for 
250Fm. For 746Cm and ?°2No, only a shallow minimum 
develops along the $32 shape degree of freedom. It was 
also shown that the evolution of the non-axial octupole 
b32 effect along the N = 150 isotonic chain is not very 
sensitive to the form of the energy density functional and 
the parameter set we used [43]. 

Both the non-axial octupole parameter $32 and the en- 
ergy gain due to the (632-distortion reach maximal val- 
ues at 748Cf in the four nuclei along the N = 150 iso- 
tonic chain [43]. This is consistent with the analysis given 
in Refs. [42,45] and the experimental observation that in 
48Cf, the 27 state is the lowest among these nuclei [41]. 
These results indicate a strong Y32-correlation in these nu- 
clei. 


4. Summary 


In this contribution we present the formalism and some ap- 
plications of the multi-dimensional constrained relativistic 
mean field (MDC-RMF) models in which all shape degrees 
of freedom (6), with even u are allowed. The potential 
energy surfaces (curves) of actinide nuclei and the effect 
of the triaxiality on the first and second fission barriers 
were investigated. It is found that besides the octupole 
deformation, the triaxiality also plays an important role 
upon the second fission barriers. The non-axial reflection- 
asymmetric 832 shape in N = 150 isotones were studied 
and rather strong non-axial octupole Y32 effects have been 
found in 248Cf and ?°°Fm which are both well deformed 
with large axial-quadrupole deformations, G29 ~ 0.3. 
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